The equation of state of solid nickel aluminide 
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The pressure- volume-temperature equation of state of the intermetalhc compound NiAl was cal- 
culated theoretically, and compared with experimental measurements. Electron ground states were 
calculated for NiAl in the CsCl structure, using ab initio pseudopotentials and density functional 
theory (DFT), and were used to predict the cold compression curve and the density of phonon states. 
It was desirable to interpolate and smooth the cold compression states; the Rose form of compression 
curve was found to reproduce the ab initio calculations well in compression but exhibited significant 
deviations in expansion. A thermodynamically-complete equation of state was constructed for NiAl, 
which overpredicted the mass density at standard temperature and pressure (STP) by 4%, fairly 
typical for predictions based on DFT A minimally-adjusted equation of state was constructed by 
tilting the cold compression energy-volume relation by ~7GPa to reproduce the observed STP mass 
density. Shock waves were induced in crystals of NiAl by the impact of laser-launched Cu flyers and 
by launching NiAl flyers into transparent windows of known properties. The TRIDENT laser was 
used to accelerate the flyers, 5 mm in diameter and 100 to 400 /im thick, to speeds between 100 and 
600 m/s. Point and line-imaging laser Doppler velocimetry was used to measure the acceleration of 
the flyer and the surface velocity history of the target. The velocity histories were used to deduce 
the stress state, and hence states on the principal Hugoniot and the flow stress. Flyers and targets 
were recovered from most experiments. The effect of elasticity and plastic flow in the sample and 
window was assessed. The ambient isotherm reproduced static compression data very well, and the 
predicted Hugoniot was consistent with shock compression data. 



PACS numbers: 62.20.-x, 62.50.+p, 64.30.-|-t 

I. INTRODUCTION 

The equation of state (EOS) relating the pressure, 
compression, and temperature of solids is important to 
understand the structure of rocky planets and the re- 
sponse of materials for dynamic loading during impacts 
or energetic events such as explosions. The EOS also 
provides a test of our knowledge of underlying physics, 
particularly the states of electrons under the influence 
of the ions, described by many-body quantum mechan- 
ics. Accurate and thermodynamically-complete EOS are 
valuable as it is extremely difficult to measure the tem- 
perature of material in a shock wave experiment, and 
temperature is a key parameter determining the rate at 
which plastic flow, phase transitions, and chemical reac- 
tions occur. 

Several approaches have been devised to calculate the 
EOS essentially from first principles. One way of classi- 
fying the approaches is into those in which the electrons 
are treated explicitly, and those in which the effect of 
the electrons is subsumed into effective interatomic po- 
tentials. Although interatomic potentials can be derived 
from calculations in which the electrons are included 
explicitljii, it can be difficult to represent the non-local 



aspects of the electron wavefunctions faithfully via an 
interatomic potential, and additional complications arise 
in deriving cross-potentials between dissimilar types of 
atom, which are necessary for calculations of compounds 
and alloys. With the electrons treated explicitly, some 
EOS calculations have been made for elements of low 
atomic number in which exchange and correlation in the 
electron wavefunctions have been treated rigorously us- 
ing quantum Monte-Carlo techniques^, but these calcu- 
lations are computationally intensive and have spanned 
relatively narrow ranges in state space. Almost uni- 
versally, exchange and correlation between the electron 
states have been treated using variants of the Kohn-Sham 
density functional theory (DFT)^iii^i^, which allows the 
electrons to be represented efficiently as single-particle 
states. Similarly, excitations of the electrons and ions 
have been represented with varying rigor. The electronic 
heat capacity may be treated with simple approxima- 
tions such as the Sommerfeld model^, by populating the 
density of electron energy levels calculated at zero tem- 
perature, with a density of levels which is calculated con- 
sistently with the excitations^"^, or it may be ignored al- 
together in many situations sampling states from room 
temperature to an electron-volt or so. Ionic motion has 



2 



been incorporated through a Griineisen model (usually 
fitted to the zero-temperature energy- volume relation)^, 
by performing simulations of the classical motion of the 
atoms under the action of interatomic potentials or forces 
from the electron states-'^^, and by calculating the nor- 
mal modes of the crystal lattice and populating them as 
phonons. 

We have previously predicted EOS and phase dia- 
grams for elements using electron ground states cal- 
culated using DFT, electron excitations into the zero- 
temperature band structure, and phonon modes calcu- 
lated from electronic restoring forces as atoms are dis- 
placed from equilibrium™. An attraction of the method 
is that compounds and alloys can in principle be treated 
in exactly the same way as elements. We have subse- 
quently demonstrated that the electron ground states can 
be found, and EOS estimated using a Griineisen treat- 
ment of the thermal excitations, for the stoichiometric 
alloy NiTi^'^. A careful study of the ab initio phonon 
modes has been made at zero pressure for the stoichio- 
metric alloy NiAli^. Here, we report a more rigorous cal- 
culation of the EOS for NiAl using ab initio phononsi^ 
and electronic excitations, and comparing with static and 
shock compression data. We also investigate the use of 
an empirical relation for the zero-temperature isotherm 
of NiAl. 

Shock wave data have been obtained mainly from the 
impact of flyers launched by the expansion of compressed 
gases. In our experiments, the flyer was accelerated by 
the expansion of a confined plasma, heated by a laser 
pulse. The interpretation of shock experiments often 
depends on the properties of other materials in the as- 
sembly, such as transparent windows, and these may be 
affected by time-dependent phenomena such as plastic 
flow. We discuss the validity of our laser-flyer experi- 
ments for measuring the EOS of NiAl. 

II. THEORETICAL EQUATION OF STATE 

A thermodynamically complete EOS can be expressed 
as any thermodynamic potential with respect to its two 
natural variables, such as specific internal energy e in 
terms of specific volume v and specific entropy s. The 
atomic properties of matter lead naturally to expressions 
for contributions to the internal energy in terms of the 
volume (or mass density p = 1/v) and temperature T. 
Following our previous study of the EOS of elementsi^, 
e was notionally split into the cold compression curve 
Cc, lattice-thermal contribution e/, and electron-thermal 
contribution Ce- 

e(p, T) = e,(p) + Clip, T) + e,{p, T). (1) 

The relation e(p, T) was then used to construct the ther- 
modynamically complete relation e(s, v) and hence other 
desired quantities such as pressure, by invoking the sec- 
ond law of thermodynamics as described later. The cold 
compression energy Cc was calculated from the ground 



state of the electrons with respect to stationary atoms, 
the lattice-thermal energy e; from the phonon modes of 
the crystal lattice, and the electron-thermal energy from 
the zero-temperature electron band structure though in 
practice it was a small contribution to the EOS in the 
solid regime. 

A. Electron ground state calculations 

Over the density range of interest, the core electrons 
are affected relatively little by each atom's environ- 
ment, so pseudopotentials were used in preference to all- 
electron calculations. The electronic ground states were 
calculated quantum mechanically using the CASTEP 
computer code. This code implements the plane wave 
pseudopotential method to solve the Kohn-Sham equa- 
tions of DFT— i^i^i^ with respect to the Schrodinger Hamil- 
tonian, to calculate energies and Hellmann-Feynman 
forces and stresses. 

Pseudopotentials were generated using the Troullier- 
Martins method^^ with Is2s2p (Al) and Is2s2p3s3p (Ni) 
electron shells treated as core. Consequently 13 valence 
electrons per unit cell are considered. NiAl was treated 
as non-magnetic. 

For a cubic material, a sequence of constant volume 
calculations suffices to determine the cold compression 
curve. The volume was held fixed while restoring forces 
were calculated to determine the dynamical matrix^^.. 

At standard temperature and pressure (STP), NiAl 
adopts the CsCl (or B2) structure, with lattice pa- 
rameter 2.8870 ± 0.0001 Ai^, giving a crystal density 
of 5.912 g/cm"^. A plane wave cutoff of 900 eV and a 
symmetry-reduced 10"^ fc-point grid^^ was sufficient to 
converge the Pulay-corrected ground states22i2ii2^i2^ to 
'--^lOmeV/A'^ or better. This level of convergence in pres- 
sure was 1-2 orders of magnitude smaller than the dis- 
crepancy from measured pressures introduced through 
the use of DFT. 



B. Isotropic compression 

Isotropic compression calculations were performed to 
determine the cold compression curve edp)- The lattice 
parameter was varied between 2.0 and 5.0 A at intervals 
of 0.1 A, with additional calculations at intervals of 0.05 A 
between 2.5 and 3.1 A and 0.01 A between 2.7 and 3.0 A. 
(Figs n and H) 

It is instructive to fit the electronic structure calcu- 
lations to the Rose functional form^'', which has been 
found empirically to describe the compression behavior 
of many elements. A functional fit such as the Rose form 
is useful in that it provides a much more compact rep- 
resentation of the cold curve than does the tabulated re- 
lation obtained from a series of electronic structure cal- 
culations. It is interesting to assess the accuracy of this 
functional form in reproducing the expansion region as 
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well as compressed states. For use with atomic struc- 
ture calculations, in which the binding energy contains 
a contribution from electrons to atoms as well as atoms 
to form the solid, the variation of binding energy with 
scaled compression can be described by 



eo - ei(l + a + O.Oba^ 



(2) 



which is a slight generalization of the original form of the 
relation-^. The parameter a is a scale length defined with 
respect to the Wigner-Seitz radius rws as a function of 
specific volume. 



a{v) 



I 



(3) 



The pressure along the cold compression curve is then 



(v/vnV/^ - 1 

p = -3^0 ^^/. (1 - 0.15a + 0.05a2)e-^ 



(4) 



The zero-pressure bulk modulus is related to the binding 
energy scale by 



eiTTla 
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(5) 



The cold curve for a material is then described by three 
principal parameters - po ~ 1/vq, ei or Bq, and I - and 
an energy offset cq. The Rose function has a theoretical 
connection with simple metals, in which the compression 
properties are dominated by the electron density. Here, 
it was regarded as a convenient functional form for inter- 
polating and smoothing data. 

Non-linear optimization was used to obtain parameters 
from the quantum-mechanical predictions of the frozen- 
ion cold curve in pressure-volume or energy-volume 
space. Some optimization schemes were unstable when 
the equilibrium volume vq was included as a free parame- 
ter. In practice, ■^o was first estimated by inspection, and 
the remaining parameters {Bq, Z, and bq) were calculated 
for fixed vq. This process was repeated for a few differ- 
ent values of vq to investigate the sensitivity, and then a 
robust but inefficient Monte-Carlo optimization scheme 
was used to locate the optimum value of vq. 

The Rose function was not able to reproduce the shape 
of the quantum mechanical cold compression curve over 
its whole range. Separate fits were produced using the 
whole set of quantum mechanical states, and using only 
the states under compression, which is the more relevant 
part of the curve for the shock wave studies of interest. 
The fit to the compression portion was good: generally 
within the numerical scatter of the quantum mechani- 
cal states. It is possible that the Rose function performs 
worse in expansion because it does not represent ade- 
quately the localization of electrons on individual atoms 
as the rarefied metal ceases to conduct. The predicted 
parameters for NiAl were bracketed by the (observed) 
parameters for Ni and Al, but not so as to suggest any 
systematic trend (e.g. averaging) that could be used to 
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FIG. 1: Quantum mechanical predictions of cold compression 
curve (energy), and fit with Rose function. 
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FIG. 2: Quantum mechanical predictions of cold compres- 
sion curve (pressure), and fit with Rose function (logarithmic 
scale) . 

The absolute value of pressure is plotted in regions of 
tension. 



predict alloy properties; this highlights the inaccuracy 
of using mixture models as is commonly attempted to 
predict the properties of alloys, particularly intermetallic 
compounds^e . (Table U and Figs [1] to H) 



C. Thermal excitation of the electrons 

At elevated temperatures, excitations of the electronic 
states can contribute to the free energy. The electron- 
thermal contribution ee(p, T) was calculated from the 
Kohn-Sham electron band structure by populating the 
calculated states according to Fermi-Dirac statistics, us- 
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TABLE I: Rose parameters fitted to ab initio frozen-ion cold compression curve for NiAl. 





Po 

(g/cm^) 


I 

(A) 


Bo 
(GPa) 


eo 

(MJ/kg) 


full set 


6.200 


0.2860 


177.8 


-1069 


compression only 


6.227 


0.2966 


189.0 


-1068 


compression only, adjusted by 7GPa 


6.007 


0.2934 


164.8 


-1070 


Ni 


8.90 


0.270 


186.0 




Al 


2.70 


0.336 


72.2 





Parameters for elemental Ni and A l ' are shown for comparison. 



ing the procedure applied previously to other elements in- 
cluding Sii^. The band structure was calculated at a tem- 
perature of K in the frozen-ion approximation: the elec- 
tron states were not re-calculated self-consistently with 
the state occupations at finite temperatures. Over the 
range of states considered, which included temperatures 
up to 10000 K, the electron excitation was rather small. 
This can be seen by considering the electron-thermal en- 
ergy, found to make only a small contribution to the EOS 
in this regime, so any error in the band energies from the 
limitations of DFT or electron-phonon coupling would be 
a small correction to this small contribution. 

The electron wavefunctions were represented at a finite 
set of positions in reciprocal space, {ki}, reduced by sym- 
metry operations so states at each /c-point had a weight 
Wi. The energy levels were used directly in estimating 
ee(w,T), rather than collecting them into a numerical 
distribution function g{E). 

Given the set of discrete levels for each compression, 
the chemical potential ^ was found as a function of tem- 
perature T by constraining the total number of valence 
electrons N: 



E 



g{E.-tJ.{T))/kBT _j_ I 



(6) 



using an iterative inversion algorithm. Once //(T) had 
been determined in this way, the expectation value of 
the electronic energy was calculated: 



{E{T)) 



EiWi 



(7) 



Repeating this calculation for each compression, and di- 
viding by the mass of the atom to obtain the specific 
energy of excitation, the electron-thermal contribution 
ee(p, T) to the EOS was obtained. 

The electron-thermal energy was at least an order of 
magnitude smaller than the lattice-thermal energy over 
the range of temperatures and compressions considered. 
At a fixed temperature, the energy decreased with com- 
pression as the bands broaden with respect to energy, 
decreasing g{^). The electron-thermal energy did how- 
ever increase more rapidly with temperature than did 
the lattice-thermal energy, so at temperatures above a 
few electron-volts the electron-thermal energy dominates. 
(Fig.H) 



D. Phonon modes 

The thermal motion of the atoms and its contribu- 
tion to the EOS were calculated in terms of the phonon 
modes of the lattice, following the method described pre- 
viously for Sii^. In general, the restoring force for dis- 
placement each atom from equilibrium may be anhar- 
monic. There is a large increase in computational com- 
plexity in calculating derivatives in the lattice potential 
energy beyond those required for harmonic phonons. In 
the present work, we have calculated effective quasihar- 
monic phonon modes for relevant amplitudes of atomic 
displacement, and assessed the sensitivity of the EOS to 
different choices of displacement. 

Quasiharmonic phonons were calculated using a force- 
constant method^' Ah initio elements in the dy- 
namical matrix were obtained by calculating the force 
on all the atoms when one atom was perturbed from its 
equilibrium position, from the charge distribution in the 
quantum mechanical description of the electron ground 
state. 

When atom i is perturbed by some finite displacement 
Ui from its equilibrium position rbi, the force on each 
atom can be found from the electron ground state. If 
$ is the electron ground state energy, the calculation 
provides elements of the stiffness matrix D, between the 
displaced atom i and all other atoms j in the supercell. 



92$ 

duiduj 



d<^{aui) 1 
dui a 



(8) 



d^{aui)/duj is the force fj on atom j when atom i is 
displaced by a distance a in the direction (unit vector) 
iii. Because partial differentiation operators commute for 
a smooth function, a row and column of the eigenprob- 
lem can be determined from an electron ground state 
calculation with a displacement Ui along one of the co- 
ordinate directions, by dividing the force on each atom 
by the displacement. Other elements can be generated 
using symmetry: for the CsCl it was sufficient to perturb 
the Ni atom at (0,0,0) and the Al atom at (i, i, i) in 
turn along the [100] direction in order to obtain the entire 
matrix of force constants. 

The squared phonon frequencies are then the eigenval- 
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(9) 



where nii is the mass of atom i. The matrix on the right- 
hand side is the dynamical matrix, D: 



[D(k)]„^ ^ ^ 



d[Ui\ad[Uj]p' 



Jk.{ri~fj)_ 



1 



(10) 



where square brackets are used to denote an element of 
a matrix. 

Restoring forces were calculated and the dynamical 
matrix constructed for several magnitudes of the dis- 
placement of the atoms. This allowed the deviation from 
a harmonic potential to be estimated, and can be used 
to determine some of the anharmonic components in the 
potential. The displacements chosen were 0, ±0.001, and 
0.01 of the lattice parameter of the 2x2x2 supercell 
used. The restoring forces were calculated for a = 3.1, 
2.9, 2.7, 2.5, 2.3, and 2.1 A for the CsCl primitive ceU. 
The potential was found to be linear in the magnitude 
of displacement, to around ±0.005 eV/A. This linearity 
indicates that the potential surface experienced by the 
atoms is harmonic, an assertion which was investigated 
by constructing EOS based on different magnitudes of 
displacement. 

To calculate the density of phonon states g{uj) for 
a structure with a given set of lattice parameters, the 
phonon eigenproblem - diagonalizing D(fc) - was solved 

for each of a set of wavevectors k. These were chosen 
randomly with a uniform distribution over the Brillouin 
zone, by choosing the components to be three indepen- 
dent random numbers with uniform distribution between 
and 1. As was found previously the density of states 
converged slowly with the number of wavevectors. We 
did not attempt to calculate fully-converged densities of 
states, as previous experience showed that, as an inte- 
grated property of the density of states, the EOS was 
only sensitive to low moments of the distribution. These 
moments did not change significantly with compression. 

The variation of lattice-thermal energy with tempera- 
ture was found by populating the phonon modes accord- 
ing to Boltzmann statistics^. At a given mass density, 
the lattice thermal energy is 
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^huJi/kT _ I 



(11) 



from which the lattice-thermal contribution ei{v,T) to 
the EOS was found by normalizing to 3 modes per atom. 
(Fig.H) 

Although the phonon modes could be calculated reli- 
ably around zero pressure as has been found previously^S, 
some negative eigenvalues of the dynamical matrix were 
found at higher compressions, giving a small proportion 
of imaginary phonon frequencies. Imaginary frequencies 
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FIG. 3: Thermal energies. The lattice-thermal energy was 
calculated from phonon densities of states which were deduced 
using different atomic displacements (solid: displacement of 
0.001; dashed; 0.01); 3.1, 2.9, 2.7, 2.5, 2.3, and 2.1 A mov- 
ing up the vertical axis as the zero-point energy increases 
with compression. The electron-thermal energies were cal- 
culated by populating the zero-temperature band structure. 
Note that the electron-thermal energy decreases with com- 
pression, as the bands move to higher energies. 

Note the finite and compression-dependent value of the 
Debye temperature, at which all curves asymptote to the 
same straight line. 



indicate possible instabilities in the lattice: the restoring 
force on a displaced atom does not increase with displace- 
ment, so the displacement can grow. This situation may 
be related to a major phase transition, a minor perturba- 
tion of the structure which is dynamically stabilized at fi- 
nite temperature^S, a restoring force that is non-linear in 
displacement, or it may be caused by minor numerical in- 
consistencies when some force components which should 
be identical by symmetry are calculated redundantly in 
different directions or by displacing different atoms. In- 
finitesimal displacements will not distinguish between 
these possibilities, though finite displacements may be 
used to investigate the structure of the potential field ex- 
perienced by each atom. In either case, the local poten- 
tial is no longer quadratic, and the quasiharmonic model 
is not strictly valid. Minor perturbations of the structure 
and non-linear restoring forces can be treated formally 
by rcnormalization of the phonon modest, though this 
procedure is unwieldy to apply consistently with the de- 
tailed shape of the restoring force at finite displacements. 
In constructing the EOS for Si^^, where nonlinear con- 
tributions to the restoring force were more pronounced, 
the principal shock Hugoniot was found to be insensitive 
to details of the phonon density of states, including the 
treatment of imaginary frequencies. We investigated the 
sensitivity of the NiAl EOS to different methods for tak- 
ing account of the imaginary frequencies: treating them 
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as freely translational - contributing a heat capacity of 
^3/2 at all temperatures - or ignoring them and renor- 
malizing the density of states to 3 modes per atom. As 
shown below, the difference in shock properties was neg- 
ligible for pressures up to 70 GPa. 



E. Construction of the equation of state 

For computational convenience, the zero-point energy 
of the lattice modes was subsumed in the lattice-thermal 
contribution rather than the cold curve; thus the cold 
curve strictly assumed frozen ions. This assumption does 
not affect the accuracy of the resulting EOS, but the 
cold curve is not physically real as defined. The electron- 
thermal energy can be predicted from band structure cal- 
culations closely related to the electron ground state cal- 
culations. However, we have found for many materials 
that the electron-thermal energy makes a negligible con- 
tribution to the EOS of the solid, and the contribution in 
NiAl has been calculated to be relatively unimportant'^^, 
so it was ignored in the present work. 



1. Thermodynamic completion 

Given the total specific energy e(T) along each iso- 
chore, the specific entropy s was found by integration of 
the second law of thermodynamics {de = Tds — pdv) : 



TABLE II: Pressure offset and bulk modulus for ab fere initio 
equations of state. 



siT) 



^ dT' de 
Y^dT' 



(12) 



The (specific) free energy / was then calculated from 
/ = e — Ts, and the pressure p was calculated by differ- 
entiating /. The thermodynamic functions were repre- 
sented by tables, so local polynomials were fitted through 
adjacent sets of points - generally quadratics through sets 
of three points - to allow differentiation and integration 
to be performed. 

EOS were generated in SESAME 301 format^^. This 
consists of rectangular tables of pressure p (GPa) and 
specific internal energy e (MJ/kg) as functions of tem- 
perature T (K) and density p (g/cm'^). For each of the 
densities in the original cold curve calculations, states 
were calculated along an isochore from T = to 10000 K. 
Although the most straightforward form to generate, 
SESAME table 301 is rather inconvenient for hydrody- 
namic calculations. A standard hydrocode requirement 
is for an EOS of the form p(p, e). To derive this from a 
table 301 requires T to be found by inverse interpolation 
given p and e, and then T and p used to find p. In the 
present work, bilinear interpolation was used in finding 
p and e at states between the ordinates of p and T. 



atom 

displacement 


offset 
(GPa) 


bulk modulus 
(GPa) 


0.001 
0.01 


-7.21 
-6.65 


165.6 
154.0 



2. Adjustment to reproduce STP mass density 

As described previously^^, any EOS can be corrected 
to reproduce the observed STP state by adding a pres- 
sure offset Apc and for consistency a corresponding en- 
ergy tilt Acc = —vApc- We have found that the pressure 
offset improves the agreement with compression data bet- 
ter than do other types of adjustment, such as scaling 
of the mass density. However, the quantum-mechanical 
cold curves asymptoted correctly toward zero pressure as 
the lattice parameter became large, so a constant pres- 
sure offset would not be accurate at low density. The 
pressure offset was calculated from the pressure given by 
each EOS at 293 K and the observed STP crystal density 
of 5.912g/cm3 (Table |TT|. A modified Rose fit was cal- 
culated for the corrected cold curve (Table Of course, 
this fit does not reproduce the STP mass density because 
thermal expansion is not included. One straightforward 
test of an EOS is to calculate the bulk modulus. For NiAl, 
the bulk modulus has been measured to be 156±3GPa2i 
and as inferred from the elastic constants is 166.0 GPa'^-. 
The adjusted first principles - ab fere initio - EOS gave 
bulk moduli which lay close to these values (Table lll|) . 



III. ISOTHERMAL COMPRESSION 

Diamond-anvil measurements have been reported for 
the isothermal compression of NiAl to 25 GPa^. Isother- 
mal compression predictions were extracted directly from 
the p{p, T) table. Isotherms from the EOS based on the 
ground state energies passed through the diamond-anvil 
data within its scatter. (Figs |4] and [5l) 



IV. SHOCK COMPRESSION 

The theoretical EOS were used to predict the prin- 
cipal shock Hugoniot, by solving the Rankine-Hugoniot 
equations''^ linking the states on either side of the shock 
to its velocity Ug'- 



2P- 
Vo — 



Po 



Up = V[{P~Po){vo -v)] 



eo + -^{p + Po){vo - v) 



(13) 
(14) 
(15) 
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FIG. 4: Comparison between calculated ambient isotherms 
and diamond anvil data. Alternative theoretical EOS gave 
isotherms which were not significantly diflFerent. 
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FIG. 5: Comparison between calculated ambient isotherms 
and diamond anvil data (detail at lower pressures). 

where subscript '0' denotes material ahead of the shock 
(with Up = 0). Given the EOS in the form p(v,e) this 
set of equations can be closed, allowing the Hugoniot 
(locus of states reached from the initial state by a single 
shock) to be calculated. The phonon and electron modes 
were populated at each mass density and temperature in 
the construction of the EOS, so the Rankine-Hugoniot 
equations were in effect solved self-consistently with the 
population of the thermal modes, though this was done 
indirectly through the use of the EOS. 

Like the EOS itself, the shock Hugoniot was calculated 
up to pressures of several hundred gigapascals. The com- 
parison with isothermal compression data demonstrates 
the accuracy of the EOS, and the electron band-structure 
remains valid until the pseudopotentials start to over- 



FIG. 6: Comparison between principal shock Hugoniot from 
different phonon treatments, in density - pressure space. The 
atom displacement is expressed with respect to the superlat- 
tice parameter. 



lap, which typically requires pressures in the terapascal 
regime. As well as providing an a priori prediction of 
the EOS and Hugoniot, the high-pressure calculations 
serve to explore the sensitivity of the theoretical predic- 
tions to different assumptions and simplifications used in 
constructing EOS, and therefore where greater care - or 
experiment - is needed to constrain the EOS. 

The quasiharmonic EOS showed some sensitivity to 
the atom displacement and the treatment of the imag- 
inary modes. The variation is an indication of the un- 
certainty in the theoretical EOS, and of regimes in which 
anharmonic terms in the lattice-thermal energy (phonon- 
phonon interactions) may contribute significantly. In 
pressure-density space, the Hugoniots were very similar 
below around 80 GPa, then differed by around 15% at 
higher pressures. In shock speed-particle speed space, 
the Hugoniots varied by 5% at low pressures, became co- 
incident for shock speeds around 7km/s (60 GPa), then 
deviated by around 5% at higher pressures. The devia- 
tion was most pronounced in pressure-temperature space: 
around 15%. With the density of phonon states corrected 
to remove imaginary modes - labeled as 'rescaled' in the 
figures - the principal Hugoniot for a displacement of 
0.001 lay much closer to that for 0.01. (Figs E] to [3) 

The EOS showed little sensitivity to the inclusion of 
the electron-thermal contribution, except in pressure- 
temperature space where it rose linearly with temper- 
ature, reaching around 10% in pressure by 10000 K. The 
exaggerated effect in temperature is not surprising: over 
the regime studied the EOS is dominated by the cold 
compression curve, so a given, relatively small, difference 
in thermal pressure equates to a considerably larger dif- 
ference in temperature. (Figs [51 to [TT]) 
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FIG. 7: Comparison between principal shock Hugoniot from 
different phonon treatments, in particle speed - shock speed 
space. The atom displacement is expressed with respect to 
the superlattice parameter. 



FIG. 9: Comparison between principal shock Hugoniot from 
different electron-thermal treatments, in density - pressure 
space. The phonons were calculated with a displacement of 
0.001 times the superlattice parameter, without reseating. 
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FIG. 8: Comparison between principal shock Hugoniot from 
different phonon treatments, in temperature - pressure space. 
The atom displacement is expressed with respect to the su- 
perlattice parameter. 



FIG. 10: Comparison between principal shock Hugoniot from 
different electron-thermal treatments, in particle speed - 
shock speed space. The phonons were calculated with a dis- 
placement of 0.001 times the superlattice parameter, without 
reseating. 



V. SHOCK WAVE EXPERIMENTS 

The experiments described here were intended to help 
validate the theoretical EOS and to cahbrate the model 
of flow stress. In many cases the samples were recovered 
after the experiments, so metallographic analysis could 
be performed in future if desired. 

To generate shock states by impact, the TRIDENT 
laser was used to launch flyer plates, each of which then 
impacted a stationary target. Two series of experiments 
were performed. In the first, the flyer was Cu and the 
target comprised a crystal of NiAl, sometimes releasing 



into a polymethyl methacrylate (PMMA) window. In the 
second series, the flyer was NiAl and the target a trans- 
parent LiF crystal. The flyer speed and surface velocity 
of the sample were measured by laser Doppler vclocime- 
try. 

The laser pulse was 600 ns long, and the flyers were 
around 50 to 400 /im thick. Velocities were a few hundred 
meters per second with a laser energy of 5 to 20 J over a 
spot 5 mm in diameter. 

The accuracy of the laser flyer technique was previ- 
ously evaluated in experiments on the EOS of Cu'^^, and 
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FIG. 11: Comparison between principal shock Hugoniot from 
different electron-thermal treatments, in temperature - pres- 
sure space. The phonons were calculated with a displacement 
of 0.001 times the superlattice parameter, without rescaling. 



has also been used to measure the EOS of NiTi^-2.. The 
experiments reported here were performed at TRIDENT 
as part of the 'Pink Flamingo' (December 2001) and 'Fly- 
ing Pig' (March 2002) campaigns. 



A. Sample preparation 

Single crystal samples of NiAl were grown from the 
melt using the optical floating zone technique. Feed-rods 
were prepared via arc melting, with excess Al included 
in the initial charge to compensate for losses during rod 
preparation and crystal growth. Crystals were grown at 
15 mm/hr along (001), starting from an oriented seed. To 
compensate for the greater evaporation rate of Al, the 
initial melt was prepared with a slight excess of Al. It is 
difficult to predict the evaporation rate precisely, and the 
crystals were slightly Al-rich. Single crystal NiAl was also 
obtained from the General Electric Corp; this material 
was closer to stoichiometry. The orientation was deter- 
mined by back-reflection Laue diffraction. Samples were 
sliced parallel to (100) and (110) planes, then ground and 
polished to the desired thickness using diamond media, 
to a 1 fim mirror finish. The grinding process imparted 
some pre-strain close to the polished surfaces. 



B. Experimental configurations 

Several different configurations of flyer experiment 
were used. In all cases, the flyer was attached to its 
transparent substrate and spaced off from the target as- 
sembly by a 'barrel' comprising a stack of plastic shims. 
The barrel was typically around 500 /im long, allowing 



target 

substrate assembly 




VISAR 
beams 



FIG. 12: Schematic cross-section of laser-launched flyer im- 
pact experiments. 



enough space for the flyer to accelerate before impact. 
(Fig. [13. 

In the initial experiments with a Cu flyer and NiAl 
target, the sample covered only half of the flyer, so the 
flyer could be seen and its velocity measured. Several 
variants were used in the design of the target assembly, 
to investigate the accuracy of data from each. In most of 
these experiments, the NiAl sample was in contact with 
a PMMA release window; in some of the PMMA win- 
dow experiments, the window was stepped to provide a 
more accurate measurement of the time at which impact 
occurred and thus the shock transit time through the 
sample. Experiments were also performed in which the 
sample was mounted on a Cu baseplate: the baseplate 
obscured the view of the flyer, but shock breakout at the 
surface of the baseplate provided a measurement of the 
pressure (from the peak free surface particle speed) and 
gave a relatively accurate measurement of the time at 
which the shock entered the sample. The baseplate de- 
sign also avoided difflculties from light reflected from the 
free surface of the PMMA windows, which sometimes 
obscured the signal from the flyer or the sample. The 
baseplate had to be relatively thin to avoid decay of the 
shock, and the sample was not recovered from experi- 
ments using this design. (Fig. [T51 ) 

The flyer speeds were low compared with the sound 
speed in the target materials, so the flnite accuracy of 
assembly - several micrometers in co-planarity or thick- 
ness - and the finite temporal registration between the 
point and line VISARs resulted in large uncertainties in 
shock transit times. The second series of experiments was 
designed to minimize the uncertainty in inferred particle 
speed and shock pressure: by using the deceleration of a 
NiAl flyer with a window of known properties (LiF), the 
Hugoniot state could be inferred from the point VISAR 
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FIG. 13: Schematic of target assembly for PMMA release and 
baseplate configurations. Center: stepped window providing 
better measurement of impact time. In all cases the flyer 
approaches from the left. 



record only, the less precise line VISAR being used as a 
coarse velocity measurement (e.g. to help count fringe 
jumps in the point VISAR record) and to verify flatness. 



C. Target assembly 




impact 
fiducial 
window 



sample 



flyer on 
substrate 



FIG. 14: Example of components of shock experiment during 
assembly (stepped release window design) . The flyer is 5 mm 
in diameter; the sample is visible through the window. 
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PMMA substrates were used, coated with layers of C, 
Al, and AI2O3, in order to absorb the laser energy and 
insulate the flyer from heating^. Cu flyers were punched 
from foils purchased from Goodfellow Corp. The foils had 
striations and machining marks which could generate in- 
terference patterns that could confuse the interpretation 
of the laser velocimetry records. The foils were polished 
manually using diamond paste to reduce the regularity of 
the marks. The material for the NiAl flyers was roughly 
triangular in shape. Rather than cutting or punching 
flyer disks before assembly, the triangles were attached 
whole to the substrates, and the drive laser punched out 
the central 5 mm to form the flyer. The material re- 
maining attached to the substrate formed a seal which 
prevented plasma from the drive escaping radially. The 
barrel was much shorter than the flyer diameter, so the 
attachment of the edges did not make the central por- 
tion of the flyer curve measurably and did not reduce its 
speed. 

Components were glued together with five-minute 
epoxy. Each flyer was clamped to its substrate while 
the glue set, to minimize the thickness of the glue layer. 
The initial viscosity of the glue was low, so the glue thick- 
ness was estimated to be negligible, essentially filling in 
surface irregularities. Components of the target assem- 
bly were clamped, and glued together by small drops at 
their corners or edges to avoid introducing any layer of 
glue between components which would change the shock 
and optical properties. (Figs [U and [TS]) 




FIG. 15: Example of assembled components (stepped release 
window design). 



D. Diagnostics 



Laser Doppler velocimetry of the 'velocity interferom- 
etry for surfaces of any reflectivity' (VISAR) type^ was 
used to measure the velocity histories of the flyer and 
the sample. A point- VISAR and a line-imaging VISAR 
were used simultaneously, the point- VISAR signal being 
recorded on digitizing oscilloscopes and the line- VISAR 
signal on an optical streak camera. The point VISAR op- 
erated at a wavelength of 532 nm with a continuous-wave 
source, and the line VISAR at a wavelength of 660 nm 
with a pulsed source ^1.5 /iS long. Timing markers were 
incorporated on the streak record, at intervals of 200 ns, 
to provide a temporal fiducial and to allow non-linearities 
in camera sweep to be removed. Relative timing of the 
point and line VISARs was deduced by comparing the 
position at which a shock wave appeared in a flyer im- 
pact experiment; the relative timing had an uncertainty 
of ~3ns. The point VISAR was generally used to ob- 
served the NiAl. 
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E. Drive beam 



The TRIDENT laser was operated in long-pulse mode, 



as in the previous flyer work- 



is^ 



using an acousto- 



optical modulator to reduce the rate at which the pulse 
intensity rose and to control its shape. The drive pulse 
was chosen to be ^600 ns long (full width, half maxi- 
mum), and was delivered at the fundamental wavelength 
of the laser: 1054 nm (infra-red). The pulses generated 
were asymmetric in time, with a long tail. 

An infra-red random-phase plate (RPP) was added to 
smooth the beam; this made a significant difference to 
the spatial uniformity. The beam optics were arranged to 
give a spot 5 mm in diameter on the substrate. The drive 
energy was quite low, so no RPP shield was included. 
The RPP collected a small amount of debris from the 
substrate, but was not significantly damaged. 
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FIG. 16: Example velocity history from a window release ex- 
periment (shot 14128). 



F. Results 

Eight experiments were performed with Cu flyers im- 
pacting NiAl, and seven with NiAl flyers impacting LiF 
windows f Table \W\ . On one shot (14138), poor fringe 
contrast meant that the reference velocity could not be 
measured by velocimetry, so it was deduced from the laser 
energy instead. 

The drive energy was measured with a calorimeter, and 
the irradiance history of the drive pulse with a photodi- 
ode. The uncertainty in energy was of the order of 1 J. 
The pulse shape was repeatable at the same and different 
energies. 

The VISAR records were used to measure the velocity 
history and flatness of each flyer. The flyers were flat to 
within the accuracy of the data over the central 3 mm, 
with a slight lag at the edges. In some experiments, we 
recorded the impact of the flyer with a window and thus 
were able to measure the flatness directly after several 
hundred microns of flight. Most of the flyers were still 
accelerating slightly at the end of the record. There was 
evidence of ringing during acceleration, but no sign of 
shock formation or spall in the flyers. 

In the release window and baseplate experiments, the 
velocity history at the surface of the NiAl generally ex- 
hibited a precursor to a particle speed of 15-30 m/s ahead 
of the main shock wave. This precursor was presumably 
an elastic wave. The rising part of the shock generally ex- 
hibited some structure, consistent with reverberations of 
the elastic wave between the surface and the approaching 
shock. The peak velocity was followed by deceleration of 
25-30 m/s (with some outliers) and reverberations, con- 
sistent with spallation and ringing. fFigs [TBI to [T51) 

In the window impact experiments, the flyer speed just 
before and just after impact were determined from the 
VISAR record (Fig.[Tll). 



1200 




1000 1100 1200 1300 1400 1500 1600 1700 1800 
time (ns) 

FIG. 17: Example spatially-resolved velocity history from a 
window release experiment (shot 14140: stepped release win- 
dow). 

The velocity from successive fringes has been displaced by 
100 m/s for clarity. 

G. Analysis 

The impact of a Cu flyer onto a NiAl target, releasing 
into a window, allowed a point on the principal Hugo- 
niot to be deduced, relying on the assumption that the 
NiAl release isentrope through the shock state gener- 
ated by the impact was equal to the principal Hugoniot 
(Fig. [20]). At the pressures of a few gigapascals gener- 
ated in these experiments, this assumption is accurate to 
a percent or better. This analysis is simplest for release 
into vacuum, where the particle speed in the sample can 
be estimated as half of the free surface speed. The dif- 
ference between this and the flyer speed is equal to the 
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TABLE III: Flyer impact experiments. 







Driver 


Target 




orientation 


thiclcness 


speed 


thickness 


speed 






(^m) 


(m/s 
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(m/s 
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release window 


















14126 


(100) 


105 


329 


± 


15 


200 


331 ± 


3 


14127 


(100) 


105 




1 

zh 


OA 


217 


390 ± 


r 




14128 


(100) 


105 


225 


± 


1 n 


217 
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baseplate 


















1/11 Qfi 


(iUUj 


250 


165 


± 


10 


398 


208 ± 


2 


14137 


(100) 


250 


200 


± 


10 


398 


260 ± 


2 


14138 


(100) 


250 


185 


± 


15 


398 


195 ± 


5 


stepped disk 


















14139 


(100) 


250 


153 


± 


2 


398 


120 ± 


10 


14140 


(100) 


250 


180 


± 


2 


398 


175 ± 


10 


window impact 


















14381 


(100) 


94 


325 


± 


3 




226 ± 


2 


14389 


(100) 


375 


102 


± 


1 




72 ± 


1 


14390 


(100) 


389 


152 


± 


3 




107 ± 


3 


14406 


(100) 


350 


141 


± 


1 




101 ± 


1 


14387 


(110) 


228 


167 


± 


2 




117 ± 


2 


14388 


(110) 


262 


167 


± 


2 




116 ± 


2 


14405 


(110) 


392 


149 


± 


2 




106 ± 


1 



The baseplates were 55 fim copper. The driver speed is the peak free surface speed of the baseplate, and the flyer speed on 
impact otherwise. The impact window was a LiF crystal, 2 mm thick, with (100) planes parallel to the impact surface. The 
'target speed' in the window impact experiments is the speed of the interface between the NiAl flyer and the LiF window 

immediately after impact. 



250 



200 



E 150 



> 



50 



1600 1800 2000 2200 2400 
time (ns) 



U3 



o 
o 

> 



350 
300 
250 
200 
150 
100 
50 



-50 



f»4t0^i,0fC/ 




200 400 600 800 10001200140016001800 
time (ns) 



FIG. 18: Example velocity history at the free surface of the FIG. 19: Example velocity histories from window impact ex- 
sample in a baseplate experiment (shot 14136). periments. 



particle speed in the flyer, and the Hugoniot of the flyer 
material can be used to find the shock pressure. Refer- 
ence Hugoniots were calculated from published EOS-^. 
Hugoniot states were obtained, with pressures between 2 
and 8 GPa. States deduced by window release and free 
surface velocity were consistent. The principal Hugoniot 
from the ab fere initio EOS passed as closely through the 



data as any unique line is likely to. (Fig. [5T]) 

A complication of the window impact experiments 
compared with the Cu flyer experiments was that the 
stress state at the interface between the sample and the 
window depends more sensitively on the elastic strain in 
both components. While the elastic constants may be 
known, the flow stress (or yield strength) is not known 
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FIG. 20: Schematic of construction used to deduce Hugoniot 
point from window release data with a flyer of known material. 
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FIG. 21: NiAl Hugoniot points deduced from the impact of 
Cu flyers, compared with the ab initio equation of state. 



for all time scales. This was not such a concern for the 
Cu flyer experiments because the flow stress of Cu and 
PMMA are much lower than for LiF'^^. Calculations were 
made using different assumptions about the flow stress, 
to assess the uncertainty and the likely behavior of LiF 
in these experiments. 

Neglecting elasticity, a measurement of the particle 
speed immediately before and after impact (uq and ui) 
can be used to determine a state on the principal Hugo- 
niot of the flyer with reference to the principal Hugoniot 




particle speed (km/s) 

FIG. 22: Hugoniot points deduced for nickel - aluminum alloy, 
compared with the ab initio equations of state. 

of the window, which must be known. Immediately after 
impact - for a time dictated mainly by the shock and 
release transit time through the flyer and window - the 
material in both components at the impact surface is at 
the same pressure and traveling at the same speed. The 
Hugoniots are expressed in pressure - particle speed {p 

- Up) space. In the shocked region, the particle velocity 
with respect to the undisturbed material in the window is 
Ml, and in the flyer uq — ui. As the same pressure exists in 
both materials, and it can be calculated from the Hugo- 
niot of the window material pw{ui), a Hugoniot state 
can be determined in the flyer: {uq — ui,pw{ui)). The 
Rankine-Hugoniot relations can then be used to calculate 
the other state parameters: mass density p, shock speed 
Us, and the change in specific internal energy e (Table HVl 
and Fig.[22D. 

More rigorously, the materials at the impact interface 
are at the same state of normal stress sn rather than 
pressure, and under uniaxial compression. By the same 
arguments as above, if the shock Hugoniot of the window 

- now considering its response as a stress tensor s, at least 
the normal component - is known, then the normal stress 
of the flyer material is thus also known. The stress com- 
ponents for a solid depend on the elastic strain and the 
elastic constants. Elastic strain can be relieved by plas- 
tic flow, which limits the deviation of normal stress from 
the hydrostatic pressure. The elastic constants of com- 
mon window materials are well-known, but the finite rate 
of plastic flow ~ which can also be considered in terms 
of a flow stress or yield stress which varies with strain 
rate - means that experiments exploring the response 
on different time scales may induce significantly different 
degrees of elastic strain for the same compression. Our 
experiments used flyers which were relatively thin com- 
pared with previous experiments, so it is possible that 
the elastic strains were somewhat greater. 

Another aspect of this complication is that we ide- 
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ally want to determine the EOS. Data from window im- 
pact experiments comprise the EOS plus the elastic shear 
stress. For a single experimental measurement in isola- 
tion, there is no way to separate the hydrostatic pressure 
(EOS) from the elastic stress, but some values can be ex- 
tracted given data from several experiments at different 
pressure, and/or knowledge of the EOS or elastic and 
plastic behavior. In the present case, relevant comple- 
mentary data include the STP elastic constants of NiAl, 
our QM EOS, and low strain-rate measurements of the 
flow stress. We considered the limiting case of no plas- 
tic flow (i.e. maximum elastic stress for a given uniaxial 
compression) , and also cases with varying flow stresses in 
the window and the NiAl flyer. 

For predicted shock Hugoniots including elasticity, the 
elastic stress was added to the Hugoniot pressure cal- 
culated from the EOS. For each state on the Hugoniot, 
the uniaxial compression i] ~ p/ po was used to calculate 
the strain tensor (using the Green-St Venant finite strain 



measure^) 



i (l/r/2 - 1) 



and hence the strain deviator 




(16) 
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The corresponding elastic stress is given by 

[c[P)Ujkl 



d[e]ki 



(17) 



(18) 



where c is the tensor of elastic constants, strictly a func- 
tion of e too. In deviatoric form. 



s = a{p, e)-p{p)I 
where p{p) is the EOS and 



= Hp)] 



ij kl • 



(19) 



(20) 



Plastic flow acts to limit the components of the elastic 
strain deviator e; in common practice it may be expressed 
as a constraint on the magnitude of the stress deviator. 

The EOS is a more general representation of the bulk 
modulus 
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Comparing the response of a crystal to uniaxial compres- 
sion along [100], with stiffness to uniaxial compression 
cii, we can define a shear modulus 
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FIG. 23: Stress states deduced for NiAl, using different as- 
sumptions about plasticity in the LiF window. 



Using the published STP elastic constants for NiAl (cn 
211.5 GPa, C12 = 143.2 GPa, C44 = 112.1 GPa)^^, ^ 
34.2 GPa. 



The shear modulus was used to calculate the devia- 
toric correction to the shock Hugoniot from the EOS, for 
uniaxial compression. These corrections were calculated 
for purely elastic response, and also for different values 
of the flow stress Y chosen to improve the match to the 
experimentally-measured shock states (Table lIVp . If the 
LiF window was treated as purely elastic then no treat- 
ment of the NiAl reproduced the experimental states. 
If the LiF was treated as elastic-plastic with flow stress 
Y = 0.36 GPa as deduced from gas gun experiments, then 
the measured states were reproduced fairly well with a 
flow stress of 0.53 GPa in the NiAl. (Figs [M] and [l5l also 
showing shock states deduced from the Cu flyer exper- 
iments. The Cu flyer data were not corrected for elas- 
tic response in the window: the correction is smaller for 
PMMA.) 

The precursor waves and the deceleration following the 
peak of the shock can be used to infer flow stresses and 
spall strengths. However, simple models parameterized 
from these data can be misleading, so the detailed con- 
stitutive behavior of NiAl will be reported separately in 
the context of microstructural response. The true (i.e. 
normal) stress to induce plastic flow at lower strain rates 
has been reported at around 1.5 GPa^, corresponding to 
a flow stress of around 0.4 GPa. One would expect the 
stresses to be at least as large on the shorter time scales 
of the laser experiments, so the flow stress inferred from 
the window impact experiments is quite plausible. 
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TABLE IV: Shock Hugoniot states deduced from window impact experiments, using different constitutive models for LiF and 
NiAl. 



shot 
shot 


particle speed (m/s) 
for different LiF models 


normal stress (GPa) 
for different NiAl models 


hydrodynamic 


Y = 0.36 GPa 


hydrodynamic 


elastic 


Y = 0.53 GPa 


(100) 

14381 
14389 
14390 
14406 


99 ± 4 
30 ± 2 
45 ± 4 
40 ± 2 


226 ± 2 
72 ± 1 
107 ± 3 
101 ± 1 


3.25 ± 0.03 
1.00 ± 0.02 
1.50 ± 0.04 
1.41 ± 0.01 


6.02 ± 0.05 
1.90 ± 0.03 
2.83 ± 0.08 
2.67 ± 0.03 


3.64 ± 0.03 
1.39 ± 0.01 
1.89 ± 0.04 
1.80 ± 0.01 


(110) 

14387 
14388 
14405 


50 ± 3 

51 ± 4 
43 ± 2 


117 ± 2 
116 ± 2 
106 ± 1 


1.63 ± 0.02 
1.62 ± 0.02 
1.48 ± 0.02 


3.09 ± 0.05 
3.07 ± 0.05 
2.80 ± 0.03 


2.03 ± 0.03 
2.01 ± 0.03 
1.87 ± 0.01 



The hydrodynamic model for NiAl was used with the hydrodynamic model for LiF. The elastic and elastic-plastic models for 

NiAl were used with the elastic-plastic model for LiF. 
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FIG. 24: Hugoniot points deduced for NiAl, compared with 
the ah initio equation of state. Points with large error bars 
are from the Cu flyer experiments. 



FIG. 25: Hugoniot points deduced for NiAl, compared with 
the ab initio equation of state. Points with large error bars 
are from the Cu flyer experiments. 



VI. CONCLUSIONS 

The frozen-ion compression curve for NiAl in the CsCl 
structure was predicted using ab initio quantum mechan- 
ical calculations of the electron band structure. The cal- 
culations reproduced the lattice parameter at p = to 
~1% with no corrections, equivalent to a discrepancy 
~ 7 GPa. The ab initio pressure- volume relation cal- 
culated using the Hellmann-Feynman theorem was not 
perfectly consistent with the ab initio energy-volume re- 
lation; this reflects the lower precision of the stress cal- 
culations. The Rose functional form, found to fit the 
energy-volume relation for compression of a wide range 
of elements, was found to fit the energy-volume rela- 
tion for compression of NiAl, but deviated significantly 
in expansion. The electron band structure was used to 
predict electron- thermal excitations. The mechanical re- 



sponse in the equation of state and shock Hugoniot did 
not alter much when the electron-thermal contribution 

was included, though the Hugoniot pressure-temperature 
relation varied by around 1% per thousand kelvin. The 
charge distribution in the electron ground states was used 
to predict ab initio phonon modes, from the forces on the 
atoms when one was displaced from equilibrium. The 
phonon density of states was quite sensitive to the mag- 
nitude of the displacement - though the restoring force 
on the displaced atom itself suggested that the effec- 
tive potential it experienced was essentially harmonic — 
but thermodynamically-complete equations of state con- 
structed using the phonon modes were not sensitive to the 
details of the density of phonon states up to ^100 GPa. 
The quasiharmonic equations of state reproduced pub- 
lished measurements of isothermal compression of NiAl 
extremely well, and were also consistent with measure- 
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ments of the shock compression. 

Laser-driven flyer experiments were performed to mea- 
sure states on the principal shock Hugoniot of NiAl, by 
impacting Cu flyers into NiAl targets and by impacting 
NiAl flyers against LiF windows. Shock transit times 
were not measured accurately enough to constrain the 
EOS. Shock states from the Cu flyer data were obtained 
with respect to the Hugoniot of Cu or PMMA. These 
states were consistent with the ab initio EOS for NiAl. 
Detailed interpretation of the NiAl flyer data depends on 
the elastic-plastic behavior of the NiAl and the LiF. Ig- 
noring these elastic contributions, the shock states were 
consistent with theoretical EOS. If the LiF was assumed 
to response elastically, the stress states deduced were im- 
plausibly high. If it was assumed to response with the 
same flow stress as observed in gas gun experiments - 
which typically explore somewhat longer time scales - 
then it was possible to reproduce the shock states by 
taking either EOS and adjusting the flow stress in the 
NiAl. Reasonable agreement was obtained for the quasi- 
harmonic EOS with a flow stress of 0.53 GPa, which is 
consistent with the range of values deduced from the 
amplitude of the elastic precursor wave in other shock 
loading experiments. The quasiharmonic EOS - which 
was found to reproduce quasistatic compression data ex- 
tremely well - is thus consistent with the shock measure- 



ments. 
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